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, Contrary to the conventional belief, it was shown that the Breit equation has the 

eigenvalues for bound states of two oppositely charged Dirac particles interacting 
through the (static) Coulomb potential. All eigenvalues reduced to those of the 
\ Schodinger case in the non-relativistic limit. 
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Abstract 



Consider the bound states of two oppositely (single) charged Dirac particles 
intercting through the (static) Coulomb potential. They can be described by 
the Breit equation [1-3] 
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where we take the specific Lorentz frame (the CM system) and the gauge 
(Coulomb gauge), aj,Pj(j = 1,2) are Dirac matrices of the particle j, while 
m and M are their masses. ty a p is the 4x4 Dirac spinor wave function (a, (5 — 
1, • • • , 4). The total energy of the system (in the CM system) is writen as 



E = ^m 2 -q 2 + \j M 2 - q 2 (2) 
where a discrete values of q coressponds to a bound state. 
It is convenient to introduce the dimensionless quantities 
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a = e 2 /hc = 1/137. 
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First, we discuss the simplest case [4]: the singlet state and equal mass case 
M = m. Let F(r) and K(r) be the radial wave functions of the large-large 
(a, (3 = 1,2) and small-small (a, (3 = 3,4) components of ^ a/3 . Then we find 
from Eq.(l) 



1 + Ap ' 
F + K = h(p)e- p/2 

and h(p) should satisfy 



(5) 
(6) 



d 2 h dh ( 2 1 
dp 2 dp\ p p(p + yp) 
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where / is the orbital angular momentum, / = 0, 1, 2, • • • . 
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The non-relativistic limit of Eq.(7) means a 2 — > there, while keeping a 2 y — ■> 
const. Ignoring terms in a 2 in Eq.(7) and noticing 1/y is of the order of a 2 , 
Eq.(7) reduces to 
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which is precisely the Schrodinger case. Therefore Eq.(7) is the proper rela- 
tivistic extension in the singlet case. 



Let us take / = 0. Near p = h(p) can be expressed as 
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h(p)=p s h(p), (9) 



where s — + — Another solution with s = — 1 — y 1 — ^ is not 
acceptable as the wave function. h(p) obeys the equation 



h" + { -1 + —— + ti 

{ P P{P + VP) J 



/i(p) can be expressed in the Taylor series near p = (|p| < ^ 



oo 

h( P ) = Y,hnP n , (11) 
n=0 

which is very much different from the Schodinger case (Laguerre polynomials) 
in the small domain < p < K Perhaps this led peoples in old days the belief 
that the Breit equation is too singular and should not be solved it directly 
(but use the perturbation to deal with it to obtain reasonable eigenvalues 
of bound states[2,3]). We should notice that h(p) is reasonably close to the 
corresponding non-relativistic solution at p ^> 1/y (see Figs. 3 and 6 shown 
below) . 

We observe that p = oo is not the regular singular point in Eq.(lO). There are 
two independent solutions near p = oo of the Eq.(lO): 



h(p)=P?{l + C j + f 2 +---} (12) 

and 

= + ^ + + (13) 

where (3 and are undetermined from Eq.(lO) near p = oo. On the other 
hand we may analytically continue the solution h(p) in Eq.(ll) valid near 
p = to p — > oo. Then we should find 

h(p) = cihi(p) + c 2 h 2 (p) near p = oo, 

where Ci,c 2 ,f3 and depend on the parameters in Eq.(lO), namely a 2 (or s) 
and y. Quantum Mechanics demands c 2 should be zero, which fixes allowed 
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values of y, and hence those of E, the energy eigenvalues. Equivalently we can 
impose the condition: After analitic continuation the solution h(p), Eq.(ll), 
near p = should smoothly join the solution fta(p), Eq.(12), near p = oo. 
This condition selects the specific discrete values of (3 (which is a function of 
a 2 and y), giving specific values of y (which must agree with the values of y 
derived from c 2 = 0). 

The Eq.(10) is known as Heun's equation, we owe to Professor Th. A. Rijken 
and his son Pieter to informed us of this information. However we have not 
yet succeeded to find the analytic form of (3 or c 2 . 

Nevertheless the perturbative consideration of Eq.(10) shows that for small a 

{3 = n- 1 -s + 0(a 4 ), 
n = l,2,..., 

where perturbation is with respect to s or a 2 , and note that 
i~0(a 2 ), 

y 

1 2 

s ~ — a . 
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To be more precise we have tried to obtain an accurate value of j3 by computer 
for the lowest two state (n = 1 and 2) in the singlet case. Changing a value 
of f3, we tried to get smooth joining of two solutions h(p) starting from p = 
and hi(p) starting from p = oo as good as possible. Table I shows our results, 
where the minimum value of 



A 



1 dh 1 dhi 



h dp h\ dp 
is found at p — p\. We see there 

(3=-s 

is reasonably good when «< 1. We found 
(3 =-s {1 + 90.5 s} for n = 1 

or 

(3 =-s {1 - 11.3 a 2 }, (n = 1) 



= (14) 
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Table 1 

The Breit equation (singlet, ground state (n = 1), equal mass M = m case) is solved 
by computer, a is the fine structure constant. However, we have computed (3 for 
various values of a. s is given by Eq. (9). is determined by computer calculation. 
A is minimized discrepancy of the logrithmic derivatives between the solution from 
p = and the solution from p — ► oo, is found at pi, see Eq. (14). 
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near a = 1/137.036. The relations between f3 versus s or a 2 are depicted in 
Figs.l, 2, 4 and 5 (for n = 1 and 2, respectively). 

We also illustrate the best fit (smallest A) of the "radial wave function" h(p) = 
(F + K) e p l 2 for the ground state of 1 S in Fig. 3 

Once f3 is known the eigenvalues E can easily be found from Eq.(lO) 



-P + | coefficient of i in Eq.(lO) j = 0, 



or 



= n-l -s + 0(a 4 ) 

where 
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Fig. 3. The ground state (n = 1) wave function h and h\ for a 2 = (137 035999765)^ anc ^ 
j3 = —s + ((5/3) = 6.65244 x 10 -6 . The horizontal axis denotes the radial parameter 
p. The vertical axis denotes values of the wave function h{p) and h\(p). The red line 
denotes h(p) from p = 0. The blue one denotes h\{p) from p = 00. h(p) and h\{p) 
are very close to the non-relativistic solution 1. 



y = and n= 1,2,3, 

aq 

The binding energy is given by 



B = 2m - E = 2m 



n=l, 2,--- 



1 + 



2(n+5/3) J 



5 (3 = (a 4 ). 
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Table 2 

The Breit equation (singlet, excited state (n = 2), equal mass M = m case) is solved 
by computer, a is the fine structure constant. However, we have computed f3 for 
various values of a. s is given by Eq. (9). (3 is determined by computer calculation. 
A is minimized discrepancy of the logrithmic derivatives between the solution from 
p = and the solution from p — ► oo, is found at pi, see Eq. (14). 

a 2 
3 
2 
1 



l 

2 

J_ 

10 

1 
100 

1 

1000 

1 

10000 

1 

(136) 2 
1 

(137F 
1 

(137.035999765) 2 

1 

(138F 
1 

20000 



— s 


(3-1 




Pi 




A 






0.5 


0.1237607 




4.4 







o 


10" 


-10 


0.2928932 


0.09529568 




4.9 







o 


io- 


-10 


0.1339746 


0.06081300 




5.1 







n 


10" 


-10 


0.06458565 


0.03739075 




5.2 









10" 


-10 


0.01257912 


0.01022230 




5.2 










io- 


-10 


0.001250782 


0.001201915 




5.2 





nu 





io- 


-10 


0.0001250078 


0.0001242386 




5.2 





rsj 





io- 


-10 


0.00001250008 


1.250836 x 10" 


-5 


5.2 





r^> 


o 


10" 


-10 


6.758241 x 10~ 6 


6.792039 x 10" 


-6 


5.2 





rs^ 


o 


10" 


-10 


6.659940 x 10~ 6 


6.694315 x 10" 


-6 


5.2 







o 


10" 


-10 


6.656441 x 10~ 6 


6.691002 x 10" 


-6 


5.2 





i-s-i 


o 


10" 


-10 


6.563769 x 10~ 6 


6.599096 x 10" 


-6 


5.2 





rs_i 


o 


io- 


-10 


6.250020 x 10~ 6 


6.287715 x 10" 


-6 


5.2 







o 


io- 


-10 



°) 



°) 



This result agrees with existing perturbative results according to M. Nio[5]. 

Next we also calculate the value of (3 for the excited state (n = 2) in the singlet 
case. The results are shown in Table II, Figs. 4, 5 and 6. 

Up to now we have discussed the simplest case, 1 S and M = m. The above 
consideration can be applied in general cases (singlet and triplet cases and 
also M 7^ m) of the Breit equation where not only the Coulomb potential but 
also the one transverse photon exchange contribution are taken into account 
[4]. Now we can conclude that the Breit equation gives eigen functions and 
eigenvalues for the Coulombian bound states, which are proper relativistic 
generalization of the (non-relativistic) Schrodinger case. 
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Fig. 4. The relatons between 0—1 versus 
— s for the singlet excited state (n = 2). 
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Fig. 5. The relatons between /3— 1 versus 
a 2 for the singlet excited state (n = 2). 




Fig. 6. The excited state 



n = 2) wave function h and hi for a 2 — ( 137 035999765)^ 
and (3 = —s + {5j3) = 1.00001. The horizontal axis denotes the radial parameter 
p. The vertical axis denotes velues of the wave function h{p). The red line denotes 
h(p) from p = 0. The blue one denotes h\(p) from p = 00. h(p) and h\{p) are very 
close to the non-relativistic solution 1 — f 
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